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Abstract 

We introduce the recently developed IR-improved DGLAP-CS theory into the showers 
in PythiaS, as this Monte Carlo event generator is in wide use at LHC. We show that, 
just as it was true in the IR-improved shower Monte Carlo Herwiri, which realizes the 
IR-improved DGLAP-CS theory in the Herwig6.5 environment, the soft limit in processes 
such as single heavy gauge boson production is now more physical in the IR-improved 
DGLAP-CS theory version of PythiaS. This opens the way to one’s getting a compari¬ 
son between the actual detector simulations for some of the LHC experiments between 
IR-improved and unimproved showers as PythiaS is used in detector simulations at LHC 
whereas Herwig6.5, the environment of the only other IR-improved DGLAP-CS QCD MC 
in the literature, Herwiril.031, is not any longer so used. Our achieving the availability 
of the IR-improved DGLAP-CS PythiaS then is an important step in the further devel¬ 
opment of the LHC precision theory program under development by the author and his 
collaborators. 



In a series of papers we and our collaborators have developed, implemented and 
applied to FNAL and LHC data the IR-improved DGLAP-CS theory in the 
Herwig6.5 environment as realized in the new Monte Carlo Herwiril.031. Because the 
IR-improvement in Herwiril.031 derives from the exact amplitude-based resummation 
theory in Refs. lOf we and our collaborators have argued [l]-|^ that Herwiril.031 should 
and does give a better £t to the FNAL and LHC data on single heavy gauge boson 
production without the need of an ad hocly hard intrinsic pt spectrum (rms pt — 2 
GeV/c) for the proton or anti-proton constituents, as one expects from observations like 
the precociousness of Bjorken scaling 13,14. As we and our collaborators continue with 


the comparisons between Herwiril.031 predictions and the recent LHCb data 15 on single 
heavy boson production and decay, we have met a matter of some concern as follows. 

In some of the spectra which depend on the transverse degrees of freedom of the heavy 
gauge bosons, detector related effects such as bin migration are based on the detector sim¬ 
ulations with the events of only some specihc MC’s and there is considerable over-head 
to re-do such simulations with Herwiril.031 events because it uses the Herwig6.5 envi¬ 
ronment whereas these detector effect modules do not use that environment currently. 
Thus, it is somewhat ill-timed to get IR-improved showers via Herwiril.031 into the LHC 
detector simulations for such effects as these important bin migration effects. We stress 
that, since the MC’s for the IR-improved and unimproved showers look very different in 
the soft regime where these migration effects tend to be more pronounced, it is impor¬ 
tant to provide a platform which will facilitate the comparison between IR-improved and 
unimproved DGLAP-CS showers in this regard. 

Accordingly, we have undertaken!^ the introduction of the IR-improved DGLAP-CS 
theory into the PythiaS environment which , at least currently, is more widely used 
in detector simulation studies at LHC. In this Letter, we describe the introduction and 
illustrate its effect on the proto-typical heavy single production px spectrum at the 


LHC. The detailed phenomenological studies will appear elsewhere 17 


Specihcally, the IR-improved DGLAP-CS theory is given in detail in Refs, [ij-j^, so 
that we will not repeat it here and we refer the reader to the latter references for its 
specihcation. We turn directly to what is needed to introduce the theory into the showers 
in PythiaS. 

Toward this end, we proceed as follows. Focusing hrst on the time-like showers in 
PythiaS, in the module TimeShower.cc we replace the usual DGLAP-CS kernels with 
the IR-improved ones in Eqs.(6) in Ref. [^. For example, whenever we have the shower 
weight factor (l-l-;^^) (note that 2 ; here is dip.z in TimeShower.cc) for a given type of QCD 
radiator color representation A with the attendant infrared point at z —>■ 1, we make the 
replacement 

(1 + Fyrs(.lA}ei‘^(l + z'Xl - zr\ (1) 

where the IR improvement exponents 'Jai^a and the YFS infrared function Fyfs(^) are 
given in Eqs.(7) and (S) in Ref. [^. 


Mhe reader interested in the chronology of the theory can see Refs. 


11 


12 


for the original Abelian 


gauge theory development and application of the approach; here, the non-Abelian generalization is needed. 
^We thank here Dr. Jesper Christiansen for useful private communications. 
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Continuing in this way, when we meet the shower weight factor (1 + z^) in the gluon, 
G, splitting to G G with the infrared point at 2 ; —)■ 1, we make the replacement 

(1 + z^) ^ FYFs(7G)e^^«(l + ^')(1 - zy^. (2) 

Finally, when we meet the shower weight factor + (1 — z^) in the splitting G ^ qq 
we make the replacement 

(z^ + (1 - zY) ^ -FYFs(7G)e!"°(z"(l - zyo + (1 - (3) 

These changes in module TimeShower.cc realize the IR-improved DGLAP-CS theory in 
time-like showers in PythiaS. 

Turning next to the space-like showers in PythiaS, we act on the module SpaceShower.ee 
as follows. When we meet the shower weight factor ( 1 —^( 1 —^))^ in the splitting G —)■ G G, 
we make the replacement 

(1 - z{l - z)y FVFs( 7 G)e 5 '^‘^ ^(1 - zfz^^ + z^{l - z^^ + - zf‘{z'^^ + (1 - ■ 

(4) 

When we meet the shower weight factor (1 -1- (1 — z^) in the splitting q —)■ G{z) q we 
make the replacement 

(1 + (1 - zf) ^ FYFs( 7 ,)e 5 ''‘' (1 + (1 - zf) z^^. (5) 

Continuing in this way, when we meet the shower weight factor (1 -|- z'^) in the splitting 
q —)■ q{z) G we make the replacement 

(1 + z^) ^ FYFs(7je^^'' (1 + z^) (1 - zy-^. (6) 

We note as well that mass corrections for the heavier quarks also receive the same IR 
improvement factors here. For the splitting G —)■ g g we make the replacement of the 
splitting weight factor {z'^ -|- (1 — 2 ;)^) as indicated above in ^ for light quarks. For 
massive quarks, the corresponding mass correction in the weight factor has its factor of 
22:(1 — z) replaced according to the rule: 

{2z{l - z)) ^ FYFs( 7 G)e^''° {z{l - z){z^^ + (1 - zy^) . (7) 


With these replacements in module SpaceShower.ee, we have introduced the IR-improved 
DGLAP-CS theory into the space-like showers of PythiaS. 

While detailed illustrations of the resulting IR-improved phenomenology will appear 
elsewhere 17 as we have noted, here we will use the pt spectrum in single heavy gauge 
boson production at the LHC to illustrate the expected size of the IR-improvement effects 
in the PythiaS environment. Accordingly, we show in Fig. [T]the px spectrum at the LHC 
for single Z/ 7 * production when the ems energy is 7 TeV and 13 TeV. We see that the IR 
improvement has the similar size effect at 7 TeV as we have seen in Herwiril.031 1^14,15 


in the Herwig6.5 environment. Here, we stress that the results in Fig. [^have no intrinsic 
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Figure 1: Comparison between IR-improved and unimproved pt spectra at the LHC 
as predicted by PythiaS for single Z/y* production at cms energies 7 TeV and 13 TeV: 
blue(green) squares correspond to IR-improved(unimproved) results for 7 TeV cms energy; 
triangles correspond to the analogous results for 13 TeV cms energy. In black and white 
print, blue(green) corresponds to dark(hght). The results presented here are untuned. 


Pt for the partons in the incoming beams. But, we note that it increases the unimproved 
PythiaS prediction in the hrst bin without the need of ad hoc manipulations as presented 
in Ref. [^. We can see from the comparisons between Pythia8(Pythia6) and ATLAS 


data in Fig. 10(Fig. 8) of Ref. 19 that the increase in the hrst bin regime is in the right 
direction to improve the agreement with the data without ad hoc parameter manipulations 
and this will be studied in more detail elsewhere 17 . These effects must be taken into 


account in analyzing the LHC data in the context of precision QCD for LHC physics, be 
it backgrounds for discoveries or SM tests. 

To sum up, we have introduced the IR-improved DGLAP-CS theory into the showers in 
PythiaS. The size of the effects are similar to those seen in Herwiril.031 in the Herwig6.5 
environment. We encourage experimentalists to use this IR-improved version of PythiaS 
to explore the interplay of IR-improvement with estimation of detector effects, especially 
when high precisions on differentially exclusive spectra are desired. The IR-improved 
version of PyhtiaS may be obtained from the website: http://bflw.web.cern.ch/ . 

In closing, we thank Prof. Ignatios Antoniadis and Prof. W. Lerche for the support 
and kind hospitality of the CERN TH Unit while part of this work was completed. We also 
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thank Profs. T. Sjostrand and P. Skands and Dr. J. Christiansen for helpfnl discnssion. 
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